Trojan states of electrons guided by Bessel beams^ 



Iwo Bialynicki-Biruljjj 

Center for Theoretical Physics, Polish Academy of Sciences 
Al. Lotnikow 32/46, 02-668 Warsaw, Poland and 
Institute of Theoretical Physics, Warsaw University 



in 
o 
o 

(N 



<Z3 

O 

-f— > 

Oh- 

q 

O 



> ■ 

in ■ 

<n : 
o . 

(N ■ 
O ' 

in : 

o . 

o : 

43 : 

Oh. 



13 



Zofia Bialynicka-Birula 

Institute of Physics, Polish Academy of Sciences, 
Al. Lotnikow 32/46, 02-668 Warsaw, Poland 



Bartosz Chmura 
College of Science, Cardinal Stefan Wyszynski University, 



Warsaw, Poland 



Previous work [I. Bialynicki-Birula, Phys. Rev. Lett. 93, 20402 (2004)] is extended to cover 
more realistic examples of electromagnetic waves, viz. the Bessel beams. It is shown that electrons 
may be guided by a Bessel beam with nonvanishing orbital angular momentum. The mechanism for 
trapping the electrons near the electromagnetic vortex line of such a wave field is the same as for the 
Trojan states of Rydberg electrons produced by a circularly polarized wave. The main difference is 
that in the present case the transverse motion of electrons in a beam is confined under the action of 
the electromagnetic wave alone, no additional attraction center is required. We also discuss briefly 
the motion of electrons in Neumann and Hankel beams. 
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I. INTRODUCTION 



The purpose of this paper is to show that Bessel beams 
of electromagnetic radiation (described in detail in the 
textbook by Stratton may serve as beam guides for 
charged particles. The confining mechanism in the trans- 
verse direction can be explained as due to an interplay 
between the Lorentz force and the Coriolis force in the 
frame rotating with the electromagnetic wave. Exact 
analytic solutions of the Lorentz, Schrodinger, Klein- 
Gordon, and Dirac equations describing beams of charged 
particles moving in the presence of an electromagnetic 
wave with a vortex line have been presented in 4] for 
a special, very simple form of electromagnetic wave car- 
rying angular momentum. These fields are not realis- 
tic since the electric and magnetic fields grow without 
bound with the distance from the vortex line. However, 
the motion of particles in such Maxwell fields helps to 
understand the confinement mechanism of particles by 
electromagnetic vortices. In addition, these simple so- 
lutions approximate very well more realistic solutions in 
the vicinity of vortex lines. In the present paper, we 
shall show that the same confining mechanism is respon- 
sible for guiding electrons inside Bessel beams of electro- 
magnetic field. Bessel beams are still not fully realistic 
because the field vectors fall off too slowly to make the 
energy finite, but they are much closer to the physical 



"This work combines two topics: Bessel beams of electromagnetic 
radiation and Trojan states of electrons. To both of these research 
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reality. 

Bessel beams of light were produced for the first time 
by Durnin, Miceli, and Eberly pj using an annular slit. 
Later jBessel beams were produced also by other methods 
HESEHIIiimill- In order to trap electrons, as 
we shall explain in the present paper, higher order Bessel 
light beams are more useful. They were produced first 
by an axicon 0] an d later in biaxial crystals 0] 



II. BESSEL BEAMS OF ELECTROMAGNETIC 
RADIATION 

Bessel beams appear in a natural way as solutions of 
Maxwell equations in cylindrical coordinates 0. These 
solutions are conveniently described using the (diff erently 
normalized) Riemann-Silberstein 0, 0, 0, ITsj vector 
F, 



F = E + icB. 



(1) 



With the use of the complex vector F, we may rewrite 
all four Maxwell equations as two equations 



id t F = cVxF, V-F = 0. 



(2) 



The separation of the complex vector F into its real (elec- 
tric) and imaginary (magnetic) parts will be needed when 
writing down the equations of motion for electrons. 

Since we are interested in the beam-like fields, we shall 
seek the solution of J5J in the form 



(3) 
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where e = ±1. Substituting this Ansatz into the Maxwell 
equations (J2J, we obtain 



(4) 



From the first two equations we may determine F x and 
F y in terms of a single complex function ip(x,y) 



F x \ 




isk z F y 




= c iek z F x 
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F, 



ckgV x ip(x, y) - iujVyipix, y) 



ck 2 ± 



ck z V y ip(x, y) + iu>V x ip(x, y) 



ck\ 



(5a) 

(5b) 
(5c) 



where k± = ^Jio 2 /c 2 — k 2 . Upon substituting these for- 
mulas into the third equation, we obtain the Hclmholtz 
equation in 2D that must be satisfied by ip 



(Vl+V 2 y + kl)Tjj(x,y)=0. 



(6) 



Every solution of this equation gives rise to a non- 
diffracting beam. Various analytic solutions may be ob- 
tained by separating the variables. 

There are three coordinate systems which allow for 
the separation of variables: polar, elliptic, and parabolic 



coordinates (cf., for example, Ref. 0])- The sepa- 
ration of variables in elliptic and parabolic coordinates 
in the Helmholtz equation leads to Mathieu and Weber 
functions, respectively. The corresponding nondiffract- 
ing beams look quite intriguing but they seem to be very 
difficult to produce in reality. In the present paper we 
shall restrict ourselves to the separation of variables in 
polar coordinates that leads to Bessel functions. In the 
degenerate case, when k± = 0, the Helmholtz equation 
reduces to the Laplace equation which separates in many 
other coordinates 0] and has a plethora of solutions. 
In particular, every analytic function of either x + iy or 
x — iy is a solution. 

The function i/j for the Bessel beam will be chosen in 
the form 



i/j(x,y) =E Q (x + iy) 



M Jm{k^p) 



P 



M 



(7) 



where Eg = cBq is the field amplitude measured in units 
of the electric field. 

The Bessel beam may be characterized by four "quan- 
tum numbers" k z ,kj_,M, and e. The meaning of these 
numbers in terms of the associated eigenvalue problems 
is discussed in the Appendix. According to Eqs. JSJ), the 
solution of Maxwell equations (0 , characterized by these 
four numbers, has the form 



{k z k ± Me} 



E e 



ie{k z z—ujt) 



/ a + (x + i y) M ~ 1 Jm -^ p) +a^(x + ty) M +^ ' h '^ p) 
ia+{x + iy) M - 1 JM -/ k t p) - ia-(x + iy) M+1 JM+ p l$t p) 
V -ie(x + iy) M 



(8) 



where u) = cy k 2 + k 2 _ 1 a± = (ui/c ± k z )/2k±, and 
p = \/ x 2 + y 2 . In the derivation, we used the formu- 
las (x\7 x + yV y ) = pd p , (xV y — y\7 x ) = 9^, and also the 
relations between the Bessel functions and their deriva- 
tives d p J M (k±p) = TMJ M (k±p)/p±kj_J MT i(k±p). 

Since Bessel beams carry angular momentum, the elec- 
tric and magnetic fields rotate as we move around the 
beam center (the z-axis). Moreover, the z-axis is at the 
same time the vortex line (except, when M = 0) accord- 
ing to the general definition proposed in 20]. 

The Bessel beam for M — 2 will play a special role in 
our analysis because it is directly related to our earlier 
work. Namely, the limit of the Bessel beam with M = 2, 
when k± — * 0, is the following solution of the Maxwell 
equations 



This electromagnetic field is a good (but not uniform) 
approximation to the Bessel beam for M = 2 in the re- 
gion where k±p and k\zjk z are much smaller than 1. It 
is the simplest example of a solution of Maxwell equa- 
tions characterized in [2lJ as "a vortex line riding atop a 
null solution" (null solution means that E 2 — c 2 B 2 = 
and E-B = 0). Electromagnetic wave (JHJ is not a plane 
wave but it has the properties found before only for plane 
waves. As has been shown in 0, one may find analytic 
solutions of the Lorentz equations of motion of a charged 
particle in this field and also analytic solutions of the 
Schrodinger, Dirac, and Klein -Gordon equations. In the 
present work we have used this exactly soluble case as 
a guide in our study of the particle's motion in a Bessel 
beam. 



(1) 



(9) 
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III. MOTION OF CHARGED PARTICLES IN A 
BESSEL BEAM 

We shall analyze the motion of a charged particle in 
a Bessel beam in a relativistic formulation in view of 
possible applications to highly energetic electrons. The 
equations of motion are in this case most conveniently 
expressed in terms of derivatives (denoted by dots) with 
respect to the proper time r = J dt^Jl — v 2 /c 2 

mX^T)=ef^(X(r))X u (r). (10) 

The trajectory is described by four functions of r 

= U(r), V (r), C(r), 9(t)} = {X(r), 9(r)}. (11) 

The equations of motion to be solved, in the three- 
dimensional notation have the form 

m \ c 
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0(t) = — »F(A>X(t). 



(12b) 



In our analysis there will always be a distinguished wave 
frequency lo and the corresponding wave-vector length 
k = bj/c. Therefore, it will be convenient to use l/u) and 
1/k as the natural units of time and distance. There are 
also the characteristic amplitudes of the electric field Eq 
and of the magnetic field Bq. Finally, we shall measure 
the velocity of electrons in units of c. We would like to 
stress that all values of electron velocities appearing in 
this paper are the derivatives with respect to the proper 
time r. They can exceed the speed of light since they dif- 
fer from the laboratory velocities by the relativistic factor 
7 = 1/yl — v 2 /c 2 > 1. In these units, the strength of 
the interaction of the electron with the electromagnetic 
field is characterized by a single dimensionless parame- 
ter a. This dimensionless parameter is known either as 
the laser- strength parameter eEo/mcu or the wiggler pa- 
rameter eBo/muj. Since in our case Eq = cBq, these 
two numbers as equal. The equations of motion for the 
dimensionless quantities have the form 



X^{r)=ar{X(r))X v {r), 



(13) 



where f^ v = f^ v /Eq is the dimensionless field and r is 
measured now in units of 1/lo (t is now effectively equal 
to tuj). In the next Section, these equations of motion 
will be solved numerically for various initial conditions 
and Bessel beam parameters. All calculations and plots 
in this work were done with Mathematica 1221. 
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IV. ELECTRONS GUIDED BY BESSEL BEAMS 

Bessel beams are capable of trapping and guiding elec- 
trons even when they have substantial initial transverse 
velocities. For example, in the optical case (27r/fc = 



FIG. 1: The trajectories of electrons guided by the Bessel 
beam (M — 2) for a = 0.0002, the longitudinal initial velocity 
0.002c, and three values of the initial transverse velocities: 
0.0004 c, 0.0008 c, and 0.0012 c. All distances are measured in 
units of c/lu. 
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FIG. 2: The trajectories of electrons guided by the Bessel 
beam (M = 2) for a = 0.14, the longitudinal initial veloc- 
ity 0.05 c, and three values of the initial transverse velocities: 
0.04 c, 0.08 c, and 0.12c. All distances are measured in units 
of c/lo. 



FIG. 3: The trajectory of a nonrelativistic electron, projected 
on the xy plane, confined by a Bessel beam. The parameters 
are the same as in Fig. 



632.8nm), studied in Ref. 0, electrons with initial trans- 
verse velocity as large as 0.0012 c are trapped by a Bessel 
beam of moderate intensity of the order of 10 14 W/m 2 . In 
Fig.^we show the electron trajectories obtained for three 
different transverse velocities and for a = 0.0002. In all 
figures presented in this paper the ratio of the transverse 
wave vector k±_ to the longitudinal component of wave 
vector k z is 1:100. 

The trapping of relativistic electrons requires higher 
values of a. This can be achieved either by increasing the 
intensity or lowering the frequency. In Fig.|21we show the 
trajectories of electrons with initial transverse velocities 
0.04 c, 0.08 c, and 0.12 c and for a = 0.14. This value of a 
may, for example, correspond to the microwave frequency 
uj = 2tt x 10 9 Hz and the intensity 1.5 x 10 9 W/m 2 (calcu- 
lated from the formula / = 3.0444 x 10~ 7 a 2 i/ 2 , where / 
is in W/m 2 and v is in Hz). In Fig.[3]and Fig.Q]we show 
the projection of the electron motion on the xy plane for 
two different sets of initial conditions. It is clearly seen 
that the motion in the transverse plane is confined to 
the vicinity of the vortex line but its details depend very 
much on the initial data. There is a substantial differ- 
ence between the slow and fast electrons — relativistic 
trajectories exhibit much more elaborate patterns. 



V. ELECTRONS TRAPPED IN HIGHER 
ORBITS 

The oscillatory behavior of Bessel functions suggests 
a possibility of trapping the electrons between two adja- 
cent maxima. Of course, a Bessel beam is far from be- 
having like a static potential. However, it does produce 
something like ring-shaped barriers in the transverse di- 
rection. To illustrate this point, we show in Fig. [S] the 
surface representing the radial component of the electric 
field. There are regions at distances of about 600, 1200, 
and 1800 units, where the electric field forms potential 
wells of a sort where the electrons can perhaps be kept 
on orbits. The calculation of the electron trajectories in 
these regions (Fig. [SJ fully confirms this expectation. We 
find there stationary (though wiggly) orbits. In contrast 
to ordinary bound states in static potentials, kinetic en- 



FIG. 4: The trajectory of a relativistic electron confined by 
the Bessel beam, projected on the xy plane. The parameters 
are the same as in Fig. |2] 




2000 

FIG. 5: The radial component of the electric field E p for the 
M = 2 Bessel beam plotted as a function of x and y evaluated 
at z = and t = 0. 



ergies of electrons trapped in Bessel beams are lower for 
higher orbits. 



VI. TRAPPING OF ELECTRONS IN 
NEUMANN BEAMS 

In addition to regular solutions, the Helmholtz equa- 
tion has also solutions with singularities. These singular 



FIG. 6: These trajectories represent the motion of an elec- 
tron with the initial velocities (— 0.5, 0.5, 0)c, (— 0.1, 0.1, 0)c, 
and (—0.15, 0.15, 0)c trapped on "higher orbits" when a — 
0.02. In order to obtain these orbits, the initial positions of 
the electron were chosen sufficiently far from the center as 
(1800,600,0), (1200,400,0), and (600,200,0), respectively. 
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FIG. 7: These trajectories represent the motion of an electron 
with the initial velocities (— 0.1, 0.1, 0)c and (—0.15, 0.15, 0)c 
trapped on "higher orbits" in the M — Neumann beam 
when a = 0.02. In order to obtain these orbits, the initial 
positions of the electron were chosen sufficiently far from the 
center as (800,400,0) and (600,200,0), respectively. 
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FIG. 8: The electric field in the Bessel beam for M — 2, 
projected on the xy plane and near the vortex line. This field 
has been evaluated at z — 0, but for all values of z it has the 
same general structure. 



solutions must be excluded if we allow the field to oc- 
cupy the whole space. However, when portions of space 
where the singularities occur, due to the presence of some 
obstacles, are inaccessible, then these singular solutions 
must, in general, be included to satisfy the boundary 
conditions. This takes place, for example, in the case of 
cylindrical coaxial lines. In order to satisfy the boundary 
conditions we have to include in the solution, in addition 
to Bessel functions Jm also the Neumann functions Ym 
(Bessel functions of the second kind). Since Neumann 
functions satisfy the same differential equation as the 
Bessel functions, the solutions of the Maxwell equations 
describing Neumann beams can be obtained directly from 
our formulas JSJ by replacing all Bessel functions by the 
corresponding Neumann functions. This will give a solu- 
tion of the Maxwell equations everywhere, except on the 
line p = 0. Assuming that the vicinity of this line is in 
some way shielded, we can study the motion of electrons 
in the region where the electromagnetic field is regular. 
In Fig. [7| we show two trajectories of electrons that were 
obtained under similar conditions as those in Fig. [6] but 
for a Neumann beam with M = 0. We can clearly see 
that the same mechanism of stabilization in the trans- 
verse plane is in place also for Neumann beams. 



VII. TROJAN MECHANISM OF ELECTRON 
TRAPPING 

In order to understand the mechanism of electron trap- 
ping near the electromagnetic vortex line, we shall rewrite 
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FIG. 9: The gradient of the saddle surface z(x, y) = x 2 — y 2 , 
projected on the xy plane. This plot is to be compared with 
the plot of the electric field in the Bessel beam Fig. |H| 



the formulas for the Bessel beams in terms of radial 
F p = (xF x +yF y )/ p and azimuthal F^ = (~yF x + xF y )/ p 
components of the vector F in cylindrical coordinates 



F, 



\ I a + J M -i(k_i_p) + a-J M +i(k±p)\ 

F =E e^\ ia+J M -i(kxp) ~ ia-J M +i(kxp) (14) 
F Z J \ -ieJ M {k±p) J 



where ip — e(k z z — cut) + M<f>. This representation ex- 
hibits clearly a screw symmetry of the Bessel beam; 
changing simultaneously z and <p in the right proportions 
z — > z + eM<po/k z and (j> — ► <f> — <fio , leaves the field F un- 
changed. It can also be shown that, as time goes by, at 
each point in space the tips of the electric and magnetic 
field vectors follow each other tracing the same ellipse 
with the frequency to. The parameters of these ellipses 
depend on p and not on z and cf>. Each ellipse lies in a 
plane determined by its normal vector N 






sJ M (kj_p) 
a+Jni-i(kj_p) - a-JM+i[k±p) 



(15) 



We can freeze the motion of the electric and magnetic 
field vectors by going to a new coordinate frame rotat- 
ing with the frequency ui/M. The relevant coordinate 
transformation in the cylindrical coordinate system has 
the form 4> — > <j> + eut/M. This transformation elimi- 
nates the time variable — in the rotating frame the elec- 
tromagnetic forces become time independent. A typical 
configuration of the electric field is shown in Fig. [8] This 
configuration resembles the forces acting on a ball mov- 
ing on the saddle surface, shown in Fig. Obviously, 
such a field configuration does not have a stable equilib- 
rium point. However, in a rotating frame, in addition 
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to the Lorentz force, there appear also the Coriolis force 
and the centrifugal force. These two inertial forces, to- 
gether with the electric field of the wave, are responsible 
for the electron trapping in the transverse direction. The 
magnetic field plays a less important role in the trapping 
mechanism, as illustrated in Figs.llOlandllll The present 
case belongs to the same category of phenomena as the 
Trojan asteroids |23|,|24|, the Trojan states of electrons in 
atoms 0, |25| or in molecules [2(| , and the Paul trap [2^ • 
In all these systems periodically changing forces lead to 
a dynamical equilibrium. In the case of Trojan states, 
the periodical changes of the forces are due to rotation. 
In the rotating frame the Coriolis force and the centrifu- 
gal force create a dynamical equilibrium in an otherwise 
unstable system. This mechanism is very well illustrated 
with the use of a mechanical model, a rotating saddle 
surface, displayed by W. Paul during his Nobel lecture 
[27j • In our case, the rotating electric field plays the cru- 
cial role. The rotating pattern of the electric field is seen 
in Fig. [H 

The same trapping mechanism operates in the case 
of Neumann beams. However, not all members of the 
family of Bessel functions can be used for trapping elec- 
trons. Special combinations of Bessel and Neumann func- 
tions — the Hankel functions 11$ = Jm + iYm and 

h'm = J M ~~ ^m — are °f interest because they describe 
outgoing and incoming waves 0- The Hankel beam is 
described by the Eq. JSJ) in which all functions Jm are 
replaced by either H M ^ or H M \ The Hankel beams do 
not seem to trap charged particles — we have not been 
able to find trapped trajectories. This is presumably due 
to a different structure of the field vectors. In Fig.lT31we 
display the radial component of the electric field for the 
M = Hankel beam. It clearly has a different charac- 
ter that in the case of a Bessel beam (cf. Fig. |3J). The 
ring-shaped barriers are now even more pronounced than 
those found for a Bessel beam shown in Fig. [SJ The lack 
of trapping, however, can be explained by a completely 
different pattern of the electric field shown in Fig. ITU 
The lines of force now spiral in, instead of forming the 
saddle pattern of Figs. |S1 and [51 



VIII. SCATTERING OF ELECTRONS OFF 
BESSEL BEAMS 




FIG. 10: The same choice of parameters as in Fig. but 
with the electric field of the wave turned off. Note that the 
amplitude of transverse oscillations is much larger than in 
Fig. EH 
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Bessel beams are the strongest at the first maximum 
and the value of the field at subsequent maxima decreases 
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field around the vortex line (cf. Fig. I15fl . 
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FIG. 12: The electric field of the Bessel beam for M = 2, 
projected on the xy plane is shown as a function of time. 
This field has been evaluated at z — for t changing from 
to 57r/12 in units of l/u>. Time is increasing when we move 
down in each column. The electric field is rotating clockwise 
at each point. 




FIG. 11: The same choice of parameters as in Fig.|5|but with 
the magnetic field of the wave turned off. The differences 
between these trajectories and the ones obtained with the 
magnetic field on (Fig.|5J is barely visible. 



FIG. 13: The radial component of the electric field E p for 
the M — Hankel beam plotted as a function of x and y 
evaluated at z — and t = 0. 
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FIG. 15: Scattering of electrons off a Bessel beam. The tra- 
jectories represent scattering of electrons in the plane perpen- 
dicular to the beam (the beam center is at the origin) with two 
initial positions and with the initial velocities (—0.02,0,0) c, 
(-0.03, 0, 0) c, (-0.04, 0, 0) c, and (-0.05, 0, 0) c. The strength 
of the Bessel beam is determined by a = 0.005. 



FIG. 14: The electric field in the Hankel beam for M = 0, 
projected on the xy plane is shown as a function of time. 
This field has been evaluated at z = for t changing from 
to 7r/2 in units of the wave period. Time is increasing when 
we move down in each column. The electric field exhibits a 
totally different pattern than the one in the Bessel beam. 
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APPENDIX A: LABELING MODE FUNCTIONS 
WITH QUANTUM NUMBERS 

The mode functions of the electromagnetic field are 
the analogues of the eigenfunctions of a set of opera- 
tors in quantum mechanics. This analogy becomes even 
more succinct when the Riemann-Silberstein vector F is 
treated as a wave function of the photon |2^. The four 
commuting operators whose eigenvalues serve as quan- 
tum numbers that characterize the mode functions are 
in our case: the z-component of the wave vector oper- 
ator (or momentum divided by the Planck's constant) 
k z = — iV z , the length squared of the transverse wave 
vector k\ = —A± = — — V^, the projection on the 
z-axis of the (dimensionless) total angular momentum 
operator J z = L z + S Zl and finally the helicity. The total 
angular momentum vector is a sum of the orbital part 
L = — tr x V and the spin part S 





-i I , ,S'„ 
i 



S z = 




(Al) 



For a plane wave the helicity is associated with the sense 
of circular polarization (right or left). More generally, 



the helicity e can be defined as the sign of the projection 
of the angular momentum J on the direction of the wave 
vector 

e = sign(J-k) = sign(S-k). (A2) 
Since the operator S k is nothing else but the curl 



(S-k) 






-v, 









-V 


-v„ 








(A3) 



we may write the Maxwell equations in the form 

id t F = c(S-k)F. (A4) 

It follows from this formula that for monochromatic 
waves e coincides with the sign of the frequency. A Bessel 
beam F{ k ^ k± M e } mav : therefore, be determined from the 
the following set of eigenvalue equations 

— iV zF{k z kj_M e} = k z F{ k ^ k±M e y, 
+ V 'y)F{k,kxM S} = k x F{k z kxMs}^ 

(-i(xV y - yV x ) + S Z J F {kzk±Me} = MF {kzk±Me} , 

sign(J-k)F {fe ^ Me} = eF {kzk±Me} . (A5) 
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